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1. Introduction 



The problem of a non-relativistic particle moving in an A^-dimensional Sombrero- 
shaped potential provides a prototype example of the spontaneous symmetry 
breaking mechanism. Yet, even when N = 1, \t is difficult to solve the corre- 
sponding Schroedinger equation with a quartic potential[l-10]. In this paper, 
we shall give explicit convergent iterative solutions for the Schroedinger equa- 
tion 

if^f = (-iv^ + V)"^ = (1.1) 

in A^-dimension and with angular momentum /. Let q be the Cartesian coor- 
dinates 

Q = {qu Q2, Qn) (1-2) 
with > 1 and the Laplacian 



The potential is 



V = ^(r^ - l)^ (1.4) 



where 

N 

r' = T.ql (1-5) 

i=l 

To illustrate our approach, it may be useful to consider first the groundstate 
(s-state) of H. By examining the shape of V, we can guess a reasonable trial 
function which approximates when is large (i.e., near r = 1). By 
taking the Laplacian of we can cast $ as the groundstate of a different 
Schroedinger equation. 
Define 

Uir)-Eo = ^-\^V^^) (1.6) 



and 



Ho^-\v'^U{r), (1.7) 



then 

Hq^ = Eq^. (1.8) 
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Introducing 

w = Ho- H = U -V (1.9) 

and 

S = Eo-E (1.10) 
the original Schroedinger equation (1.1) can be written as 

{Hq-Eq)'^ = {w-S)"^. (1.11) 

In (1.6), only the difference U — Eq is defined. The constant Eq nnay be chosen 
by requiring 

w;(oo) = 0. (1.12) 

Following Refs. [11-14], the original Schroedinger equation (1.1) will be 
solved through an iterative sequence 

[Ho - Eo)^(m) = iw- ej^im - 1) (1.13) 

with m = 1, 2, • • • and when m = 

*(0) = $. (1.14) 

The following sinnple discussion provides the motivation of this approach. Mul- 
tiplying (1.13) by $ and (1.8) by ^'(m), we find their difference given by the 
fanniliar Wronskian-type expression 

_lv • (W^'(m) - ^(m)V$) = {w- Smj^'^im - 1). (1.15) 
Integrating (1.15) over all space, we have 

em = \ wm{m - l)d^q I J W(m - l)d^q . (1.16) 
Note that if ^'^(m) is a solution of (1.13), so is 

^(m) = ^'°(m) + c (1.17) 

where c is a constant. Since $ is the groundstate of Hq, it is positive ev- 
erywhere. By adjusting the constant c, '^{m) can also be made positive ev- 
erywhere. Thus, if w is bounded, so is Em- This approach prevents the kind 
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of divergence encountered by the usual perturbative series, and enables us to 
derive explicit convergent iterative solutions, as we shall see. 

In Section 2, we show that when y is a radial potential, the successive iter- 
ative solutions ^f(m) can be solved by sinnple quadratures. For the Sombrero- 
shaped potential (1.4) the low-lying wave function is known reasonably well 
when it is in the asynnptotic region when r is large, and also when its ampli- 
tude is large (i.e., r near 1). In Section 3, a reasonable trial function $ is 
constructed whose fractional deviation from is large only for r small (i.e., 
when the absolute magnitude of is also very small); the same should also 
apply to w{r), the difference between Hq and H, given by (1.9). As will be 
proved in Section 4, w{r) is positive everywhere and its radial derivative w'{r) 
always negative, making w maximal at r = 0. These conditions enable us to 
apply the Hierarchy Theorem [12-14] which guarantees the convergence, as will 
be discussed in Section 5. In Appendix B, we show that the rate of convergence 
for large g is similar to a power series in g~^. 
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2. The Angular Momentum Operator and the Radial 
Equation 

Express the Cartesian coordinates qi, ^2, ■ ■ ■ , Qn of (1.2) in terms of the 
radial variable r and the {N — 1) angular variables 



through 



and 



with 



Ou ^2,---, On-2 and On-i (2.1) 



qi = rcos^i, ^2 = ''^sin^i cos^2, 

^3 = r sin 6*1 sin 6*2 cos 6*3, ••• , (2.2) 

qN-i = r sin6'i sin^2 • • • sin^7v-2 cos^iv-i 



qn = r sin^i sin^2 • • • sin^7v-2 sin^7v-i 



< < TT for i = 1, 2, • • • , AT - 2 
and (2.3) 
< On-i < 27r. 

Correspondingly, the line elennents are 

dr, rdOi, rsm0id02, r sin^i sin^2C^^3, 
rsin^isin^2sin^3C?^4, r sin^i sin^2 • • • sin^iV-2C^^7V-i (2.4) 



and the Laplacian operator is 

= ^ iir'^'i) - ^'^'i^ - 1) (2-5) 

where 
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sm^^ ^OidOi d9i sm^^i 

1 d , . M ^ d . 1 



C\N-2) = --^^^^^^ {sin"-' 0,—) + -^C\N - 3) 
sm O2062 062 sm 62 



(2.6) 

1 d . . ^ d . 1 



sm C7jV-2 C/"N-2 sm UN -2 



and 



£^(1) = ^ 



The square of the angular momentum operator on an n-sphere is C'^{n). From 
( 12 .61) . one sees readily that the commutator between any two £^(n) and C'^{m) 
is zero; i.e., 

[C^{n), C^{m)] =0. (2.7) 

As we shall show in Appendix A, the eigenvalues of each C'^{n) are 

/(/ + n-l) (2.8) 

with / = 0, 1, 2, Thus, for a radially symmetric potential, the wave 
function can be written as 



^(r, ^1, ^2, • • • , ^^-1) = 7^(r)e(^l, ^2, • • • , ^^-1) (2.9) 



with 



and 



c^{N-i)e = /i(/i + iv-2)e, 
jC\N-2)e = i2{i2 + N-3)e, 

c\2)e = iN-2{iN^2 + i)e 



£2(1)6 = i%_,e. 



(2.10) 
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Correspondingly, by using (2.5) we find the radial part of the Schroedinger 
equation (1.1) for the lowest eigenstate of angular momentum / to be 

+ 2^'*' + - 2) + V(t) - E ] n{r) = (2.11) 



V? = (2.12) 



with 

1 d , r)T^ d 
(r — 

J.2K ^j, ^j. ' 

and / = l\ given by the first equation of ( 12.101 ): i.e., 

C}{N - 1)6 = /(/ + iV - 2)6. (2.13) 

Define 

7^(^) =r-^?/^(r). (2.14) 

Eq . ( 12 . 1 11) becomes 
where 

k = l^K = l^\N -\) (2.16) 

and therefore k{k-\) = K{K-l) + l{l+N -2), on account K = \{N-l). 
As r ^ 0, 7?. oc and therefore il) oc r'+^ = r^. 

The radial trial function 0(r) shall be constructed in the next section. Sim- 
ilar to (2.15) it satisfies 

{Hq - Eo)(P{r) = (2.17) 

with 

The iterative equation (1.11) becomes 

{Ho - ^o)^m(r) = {w{r) - Sm)i^m-i{r) (2.19) 

where, as before, 

w{r) = U{r) - V{r). (2.20) 



Define 

fmir)=i^m{r)/(l){r). (2.21) 
As in (1.15)-(1.16), we have 

-^|:(0'^) = («^-U0Vm-i (2.22) 

and 

oo oo 

^m = J w(P^fm-idr/ J (P^fm-idr. (2.23) 



Eq. (2.22) can be readily integrated 

00 7 oo 



fm{r) = fmioo) - 2 / ^ / iw{x) - cl>\x) fm-i{x) dx (2.24) 
or, equivalently 

fmir) = fmiO) - 2 / / (Wix) - em)(p\x)fm-lix)dx. (2.25) 

Thus, with 

Mr) = 1, 

each fm{r) can be derived by quadratures. At each iteration, there is an ar- 
bitrary constant fm{oo) or /^(O) (similar to the constant c in (1.17)). As we 
shall see, different choices of this constant can lead to different convergent 
iterative solutions. 



dy ' 



Remarks: We note that in (2.11), the centrifugal potential -^1(1 + N — 2) be- 
comes the familiar P/2r'^ in two-dimension and /(/ + l)/2r^ in three-dimension. 
Correspondingly, the term k{k — l)/2r^ in (2.15) becomes (/^ — |)/2r^ in two- 
dimension, but remains 1(1 -\- l)/2r^ in three-dimension. 
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3. The Trial Function 



To construct a reasonable trial function 0(r), we start with the following two 
functions 0+(r) and (j)-{r). 

- (3.1) 



and 



where 



r + 1 r -\- a' 



0_(r) = ^(l±^)fc e-|^+^^o(r) (3 2) 

r + 1 r + a 



5o(r) = i(r-l)2(r + 2) (3.3) 



and, as before, k is related to the angular monnentum I and the dimensionality 
N by 

/c = / + i(iV-l). (3.4) 

The paranneter a is positive but free within a range, which will be specified 
below. The trial function 0(r) is given by 

_ / + 9-(l>-{r) for < r < 1 

^^^"l(p+ + P-e-iWr) for r > 1, ^^"^^ 

where Qj^ and are constants given by 

9± = 9^{-+l). (3.6) 
Thus, by construction is regular at r = 0, with its derivative 

(r-^-cj))' = ^(r"V) = at r = 0. (3.7) 
dr 

Throughout, ' denotes d/dr. By construction, 0(r) and 0'(r) are continuous 
everywhere, with r varying fronn to oo. 
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By differentiation, (J)+{r) and 0(r) satisfy respectively 

+ i^^^ -l) + V + h)c^^ = Eo0+ (3.8) 
and, as in (2.17)-(2.20), 

^~\jr^ + 2^^^^ - 1) + V + = £;o0, (3.9) 

where 

EQ = g{l + ka), (3.10) 

, kaq^ ko?q 1 /c(/c + 1) /ca . x 

^0 = , ^ + + 7 + ^ + 7 77 TT 3.11 

r(r + a) r + a (r + 1)^ 2(r + a)2 (r + a)(r + l)^ ^ 



and 



with 



w{r) = h{r) + g{r) (3.12) 



fca(l-r^)l 



where 



g(r) = j ''"^^^^r" for < r < 1 (3.13) 

for r > 1, 



3 

A = 2r/(r-y) (3.14) 



and by construction [g-/{g+e^ + 5'-)] = 9-4'-/ 4'- 
Remarks 



(i) The exponent ^gSo{r) in ( 13.11 ) and ( 13.21) satisfies 



\\±aSi(r)? = \gHr' - I f = V(r). (3.15) 

(ii) For / = and = 1, we have k = 0. Correspondingly, 0+(r) of ( 13.11 ) 
becomes 

Mr) = ^ e-^^"W (3.16) 
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which satisfies 

. 1 £ 



with 



+ y + /io)0o = (3.17) 



hoir) = j^^. (3.18) 

These equations reduce to the same expressions used before for the corre- 
sponding one-dimensional trial function[12, 13]. 

(iii) From (3.11) and (3.13), we see that at r = both h and g contain an 
r^^ pole term. This is because neither the derivative of r^^(f)+(r) nor that of 
r^^0_(r) is zero at r = 0. The trial function 0(r) does satisfy [r^^cj))' = 
at r = 0, in accordance with (3.5) and (3.7). Thus, the potential function 
w{r) = h{r) + g(r) given by (3.12) is regular at r = 0. 

Throughout the paper, we assume 

9>- + h (3.19) 

so that > 0. Correspondingly, as will be shown, w{r) is positive and has a 
discontinuity at r = 1. 
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4. Properties of w{r] 

In this section, we shall establish 

w{r) > 

and 

w'{r) < for r > 0. 
It is convenient to write ( 13 .111) as 



with 



hi(r) = 



kag. 



r[r + a) 

k{k + l) 



i=l 



h2{r) = — ■ — , hs{r) = 
r + a 



1 



hAir) = 



and h^(r) = 



(r + l)2' 
ka 



2(r + a)2 

Likewise, we decompose ( 13.131 ) for r < 1 as 



[r + a)(r + 1) 



with 



and 



9&{r) = ^— — , g7{r^ 

g+e^ + g- 

t) = 



2gg^kar 



{9+e^ + 9-){r + a] 
2kg g -a 



{g+e^ + g-)r{r + a) 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



In order to show that w = h + g satisfies ( 14.11 ). we note that for r > 1,^ = 
in accordance with ( 13. 3j) . and therefore 

w{r) = h{r). 

Since each of the hi{r) in ( 14.31) satisfies 

h^{r) > and h[{r) < 0, 
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therefore 



w(r) > and w'(r) < for r > 1. 



(4.6) 



For r < 1, it is convenient to combine first some of the hi with gj given by 
(lO) and (I431I . 
Defining 



we find 



ka 



Wj 



r[r + 
Likewise, write 



(9- 



wi = hi +§8, 



^99- 



^2 + ^7 = 9ka 
= gka 
= gka 



^J ^ kag+g-{e^ - 1) 

g+e^ + g- r{r + a){g+e^ + g-] 

2g_r 



a 



+ 



r + a ((/+e^ + (/_)(r + a) 
2(/_[(r + a) — a] 



a 



+ 



r + a (^f+e^ + 5'_)(r + a 



a 



■[1 



^9- 



+ 



2^- 



r + a' 9+^^ + 9- 9+^^ + 9- 
Separate the two terms inside the curly brackets and write 

h2+g7 = Wii + Wo 

with 

\9+e'' - 9- 



Wii = 



gko? ^" 



r + a 



g+e^ + g- 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



in which the factor inside the square brackets is identical to the corresponding 
one in ( 14. 91) . and 



Wo 



2gg^ka 

g+e^ + 9- 



(4.12) 



Next, combine the above wo with ge of (14.51) into a single term as follows 

2gg^{ka + 1) 



wvi = wo + gQ 



g+e^ + g- 



(4.13) 



13 



We now express i(; = /i + ^,forr<l,in terms of a new sum of six terms: 

VI 

w=Y.Wm (4.14) 

m—I 

with wi, wn and wyi given by (14. 91) . (14.111) and (14.131) respectively; the rest 
Win = hs, wiv = /i4 and wy = h^- In explicit forms 

Zka Ygko? 
wi = — ■ — r, wu = 



where 



r(r + a) ' r + a ' 

1 /c(/c+l) ,,,,, 

^^^^ = / ^n2 ' = ^ ^4-^^) 

(r + 1)"^ 2(r + ay 

wv = 7 — ■ — , 3nd = X25f(A;a + 1) 
(r + a)(r + 1) 



X= I- , (4.16) 



y = ^±4^. (4.17) 

^+e^ + g- 



and 



Z = J, • (4-18) 



g+g-{e^ - 1 ) 

For the with m = ///, IV and it is clear that 

Wm > and w'^ < 0. (4.19) 
For m = VI, since in accordance with (13.141) . for r < 1, A > and 

A' = 2g{l - r^) > 0, (4.20) 



wvi also satisfies (14.191 ). 
Introduce 



^ = = ^(^ - y) (4.21) 
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and 

9- 

Recall 

9+ +9- = 2^, 

on account of ( 13 .61 ). we can express g± in terms of g and b: 

ge^ , ge~^ 

9+ = TT and g. = 



cosh b cosh b 

Thus, we can rewrite ( 14.181 ) as 

Z = ^ sinh^ 

cosh 6 cosh(^ + b) ' 

and correspondingly 

Z 9 /^\/tanh^\/ cosh^ 



r cosh6Vr/V ^ /Vcosh(^ + 6)/ 

We observe that 



and for r < 1, 



^ < 

(ir Vr/ ^ 3 



|^.(l-.^)>0. 



Since for ^ positive 



, ' , (2^-smh2O<0 



and 

c?^ Vcosh(^ + b) 



4? 1^ ( — w^!^^ ) = tanh^ - tanh(^ + 6) < 0, 
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we have for r < 1 



d /Z 



dr^r 



< 0. 



Thus from ( ]4~T5l) . 



wi > 



and 



w'j <0 for r < 1. 
Lastly, we examine wu. From ( 14 .17^ we see that 

y > 0, 

and therefore, in accordance with ( 14.151) 

Ygka^ 



Wu = 



> 0. 



r + a 

By using ( 14.171 ) and ( 14.201) . we also have 

1 1 



(Inn = A^)^+e^2, 

9+& - 9- +9-' 



g2(A+26) _ I 

2ff(l-r^) 
sinh(A + 26) 



> 



for r < 1. Therefore 



{In wn)' = (InY)' 



r + a 
2g{l-r') 1 



sinh(A + 26) r + a 
1 



(r + a) sinh(A + 26) 



with 



S = sinh(A + 26) 
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and 



T = 2g{l-r'^){r + a). (4.38) 

We seek to show that T - 5 < for < r < 1. On account of (]42D|), 

{T-Sy = -4£/r(r + a) + 2£/(l-r2)-A'cosh(A + 26) 
= -4gr{r + a) + 2g{l - r^)[l - cosh(A + 26)]. 

At any r < 1, (T - 5)' < and therefore 

r(r) - S{r) < r(0) - 5(0). (4.39) 
When r = 0, A(0) = 0, 



S{0) = sinh26 



^{9+ 9-\ _ 9%- 9^- 



2^5^- 9+^ '^9+9- 
lr(^ + ! + l)^-(^-! 



1)^ 



9' 



9' 



and 



r(0) = 2ga. 

Setting the parameter a to be within an upper bound defined by 



a < 



+ 1 



9' 



fc+l)2' 

n. ' 



we have r(0) — 5'(0) < 0, and therefore 

w'u{^) < for r < 1. 
The inequality ( 14 .42^ can also be written as 



.k 1 



g- < (- + !)(- + - + !). 
a a a 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



(4.44) 
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Combining with the inequality ( 13. 191) . we require 



1<^< 



1 + -r^. (4.45) 



a 



This is the sufficient condition for w'jj < 0, and therefore also 

w' <0 for r < 1. (4.46) 
At r = 1, in accordance with (I4l1)-(l475ll 



^7(1) + ^8(1) = (4.47) 

and therefore 

gii) = m= ^.f: >0- (4.48) 

Thus, w{r) has a discontinuity at r = 1, with 

w{l-) - w{l+) = g{l) > (4.49) 

Together with ( 14 .61 ) and ( 14.461 ) . this proves that for the parameter a within the 
range ( 14451) . 

w{r) > and w'{r) < 
over the entire range of r. 

Remarks 



(i) The first inequality in ( 14.451) 



- + !<(/ 
a 

prevents the limit a ^ 0; otherwise h{r) would have a double pole r^, in 
accordance with (4.2)-(4.3). 

(ii) Choose the parameter a within the range ( 14.451) . As we will discuss in 



the next section, because of (14. Ij) . the application of the Hierarchy Theorem 
leads to a convergent iterative solution in terms of quadratures for the lowest 
eigenstate Schroedinger wave function (1.1)-(1.5) in any dimension and 
with arbitrary angular momentum I. 
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5. The Convergent Iterative Solution 



5.1 The Hierarchy Theorem 

We begin with the iterative equations (2.22)-(2.25). Recalling that /^(r) 
and 8m are the m*^ order solutions for / = i/j/cf) and 8 = Eq — E, we derive 
the corresponding m*^ order solutions for if; and E to be 

ipm = (l>fm and Em = Eq- Sm (5.1) 

Throughout this section, we assume the parameter a to be restricted by the two 
inequalities given by (4.45), Therefore, for large g, a is small. We distinguish 
two different conditions: 

{A) /m(oo) = 1 for all m (5.2) 

or 

{B) fm{0) = 1 for all m. (5.3) 
Correspondingly in (A), (2.24) becomes 




and in (B), (2.25) becomes 




In both cases, we assume that the two inequalities of (4.45) hold. Thus w{r) 
is positive and w'{r) negative at all r, in accordance with (4.1). In case (^4), 
apart from these two inequalities of (4.45) there is no other restriction on the 
magnitude of w{r). In case (B), we assume w{r) to be not too large so that 
for all m 

fmir) > at all r. (5.6) 
As we shall see, this is equivalent to requiring for all m 

/m(oo) > 0. (5.7) 
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Hierarchy Theorem[12-14] 

{A) With the boundary condition /^(oo) = 1, we have for all m 

(5.8) 

Thus, the sequences {Sm} and {/m(^)} are all monotonic, with 

ei<e2<^z<--- (5.10) 

and 

l</i(r)</2(r)</3(r)<--- (5.11) 

at all finite r. 

{B) With the boundary condition /m(0) = 1, we have for all odd m = 2n+l 
an ascending sequence 

Si<e^<e^<---, (5.12) 

but for all even m = 2n, a descending sequence 

^2>^4>^6>---- (5.13) 

In addition, between any even m = 2n and any odd m = 2/ + 1 

^2n > W- (5.14) 

Likewise, at any r, for any even m = 2n 
whereas for any odd m = 2/ + 1 

/ /2/+2(r) ' 
c^r V/2z+i(r) 

Furthermore, 

lim Em = E (5.17) 



> 0. (5.16) 
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and 

lirn^ fm{r) = fir) = ^{r)/ct>{r). (5.18) 

Thus, the boundary condition fmioo) = 1 yields a sequence, in accordance 
with (5.10), 

Ei> E2> E3> ■■■ > E, (5.19) 

with each nnennber E^ an upper bound of E, similar to the usual variational 
method. 

On the other hand, with the boundary condition fm{0) = 1, while the 
sequence of its odd members m = 2/ + 1 yields a similar one, like (5.19), with 

Ei> Es> E5> ■■■ > E, (5.20) 

its even members m = 2n satisfy 

E2<E^<Ee<---<E. (5.21) 

It is unusual to have an iterative sequence of lower bounds of the eigenvalue 
E. Together, these sequences may be quite efficient to pinpoint the limiting 

E. 

In Appendix B, by examining a simple prototype example in this class of 
problems, we shall show that for the Sombrero-shaped potential when g is 
large, the rate of convergence is similar to a power series in g~^. 



5.2 Numerical Results 



From (4.45), we see that for a given pair (/c, g), in order to apply the 
Hierarchy Theorem, the parameter a should be within a range 

0"min < < 0"max (5.22) 

with the limits amin and amax determined by 

p = 1 + (5.23) 



and 



/ = (l + -^)(l + ^). (5.24) 
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Alternatively, for a given pair {k, a), g should be within 



9min 9 ^ 9max 



(5.25) 



with 



9min 




(5.26) 



and 




k + 1 



a 



)• 



(5.27) 



Examples of these limiting values are given in Table 1. 

Figure 1 gives examples of the iterative solutions ipn 3rid En with different 
n for the parameters g = 3, k = 2, and a = 1.2. One sees that for both 
boundary conditions /n(0) = 1 and /n(oo) = 1, the convergence sets in 
rapidly in accordance with the Hierarchy Theorem. In the case of the boundary 
condition /(O) = 1, there is a limit to the range of these parameters, in order 
that /n(oo) > 0, in accordance with (5.7). For example for g = 3, the 
boundary condition /(O) = 1 can be applied only for k < 2.5, while the 
boundary condition /(oo) = 1 can be applied to any values of k. 

In Figures 2 and 3 we give the final radial wave function lZ{r) and energy 
E that satisfy (2.11) for 



Recalling (2.14) and (2.16), we have the curves in Fig. 2 for / = and 
= 3, 4, 5, and 6 corresponding to k = K = = 1, 1.5, 2 and 2.5; 
T^{r) = ■^il){r). In Fig. 3 the curves for = 3 and / = 0, 1, 2, and 3 corre- 
spond to = 1 and A; = + / = 1, 2, 3, and 4; 7^(r) = :^'^{r) = ^^'(O- 



^ = 3, / = and iV = 3, 4, 5, 6 



and also for 



(5.28) 



^ = 3, iV = 3 and / = 0, 1, 2, 3. 
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Appendix A 

Denote the solution 9 of (2.10) as 

e> = ©Lt./L(^i, ^2,---, ^A^-i). (Ai) 

In this Appendix, we shall derive its explicit form inductively. 
Write 

®M2-V;v-i(^l' ^2, • • • , ^N-i) = Zl^^s'i^i^..%_^{62, 6'3, • • • , On-i) iA.2) 

with Zi^^ depending only on 6i. By using (2.6) and (2.10), we see that Zi^^ 
satisfies 



. N-2n y^^^ ^1^) + — ::T:x7i ''W + ^ - 2j 



sin''-'^i5^i dOi sin'^i 

(^.3) 

For N > 2, '\t is convenient to denote 

2; = COS^l, ll=l, 12=171 (^-4) 

and express (A. 3) as 



2 



(l-z^)^-(iV-l).^+^(^+^_2)-!^<!I^±^ Z-^(z) = 0. (A5) 



When N = 2, the operator ^^(tv - 1) = £^(1) in the last equation of (2.6) 
is given by 

£'(1) = ^ (-4.6) 

with 6 denoting the corresponding 9n-i- Likewise, the last equation of (3.10) 
can be written as 

~e = {A.7) 

with P denoting the corresponding its solutions will be designated as 

e = zj;^ = cosW {A.8) 
for the functions even in 6, and 

e =^fj = sinW {A.9) 
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for the functions odd in 6, with 

/ = 0, 1, 2,--- . (AlO) 

For N = 3, the eigenfunction of (A. 5) are the Legendre polynomial when 
m = 0; i.e., 

zg>(.) = p,{z) = - zy {A.n) 

where / = 0, 1, 2, • • •, as before. The corresponding Zi^(z) for m > is 
given by the associated Legendre function 

zi%{z) = Priz) = (1 - z')'^j^PM (A.n) 

with 

m<l. (A13) 

In order to derive the functions Zi^{z) for AT > 3, it is useful to establish 
the following properties: 

(i) Differentiate m times the following equation for u{z) 

where a and h are constants, and denote 

v{z) = ^uiz). (A.15) 



We obtain 



with 



a^n = a + 2m 

and {A.n) 

Pm = b — ma — m{m — 1). 
(ii) Instead of (A. 14), u{z) now satisfies 

9, d'^u du , cu „ / . . 
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where a, b and c are again all constants. Write 

w{z) = {l-z^)^u{z). (A19) 



we find 



with 



(l-/)^-^„.- + B„»-^ = (A20) 



An = a — 2n 
Bn = b + n{a — n — 1) 
and (A21) 

Cn = c -\- n{a — n — 2); 

therefore, 

Bn-Cn = b-c + n. {A.22) 

We note that from (A. 11) the Legendre polynomial Pi{z) = Zi^q(z) satisfies 
(A. 14) with 

a = 2, and 6 = /(/ + l). (^.23) 
Thus, for N = odd = 2A; + 1, we have 

It can also be readily verified that for m > 0, 

Zl"j(z) = (l-z')?j^zi;\z). (A25) 
For N = even = 2k, we start from zj'^Q{z) = cosW with z = cos6, and write 

j/c— 1 

Zi"\z) = _zgi_i,o(^); (^.26) 
for m> 0, the corresponding Zi^{z) is given by the same (A. 25). 
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Appendix B 



Because of (4.1), w; > and w' < 0, and the hierarchy theorenn, the iterative 
solution with the boundary condition /(oo) = 1 is convergent for any g > 0. 
By examining a simple prototype example in this class of problems, we shall 
show that for g large, the rate of convergence is similar to a power series in 



Consider the Schroedinger equation 

1 



-^^" + iVix)-E)^ = {B.l) 



in one space dimension with 



Vix) = ?^(x' - \f (B.2) 



and 



X > 0. (5.3) 

Throughout this appendix, ' denotes ^ To simplify the analysis, we impose 
the boundary conditions 

V'(oo) = 

and at the origin (^-4) 



ip / x=0 

The trial function for the groundstate wave function is chosen to be 

(f)(x) = e-^^° (5.5) 



with 



5o(x) = i(x-l)2(x + 2); (5.6) 



it satisfies 

-^ct>" + {U{x)-g)cf> = (5.7) 

with 

U{x) = V{x) + u{x), (5.8) 
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and the same boundary conditions 0(oo) = and 



Rewrite (B.l) as 
with 

Introducing 

we find 

at X = 



5-j =^-1. (5.10) 



+ {U{x) - g)^ = {u{x) - E)^ (5.11) 



E = g-e. (B.12) 



/^t (5.13) 

-\{ct?n' = {u{x)-S)ct?f- (R14) 



/(0) = 0. (5.15) 
To fix the relative nornnalization factor between ip and 0, we innpose at a: = oo, 

/(oo) = 1. (5.16) 

The groundstate wave function il){x) of the Schroedinger equation will be 
solved by introducing the iterative sequences 

V^l(x), V^2(a:), • • • , V^n(a^), • • • 

and (5.17) 

Ell ^2-1 ' ' ' 1 Em 

with 

-\< + {U{x) - g)'^n = {u{x) - En)^n-1 (5.18) 

and 

^0 = 0. (5.19) 
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Define 

W) = %^ . (S.20) 
From the above equations, it follows that 

= {u- £:O0Vn-i; (5.21) 

furthernnore, similar to (B.15)-(B.16), 

/;(0) = {B.22) 

and 

/n(oo) = 1. (5.23) 

From (B.21)-(B.22), we have 

oo 

/ (uix) - 8n)(l)\x)fn-l{x)dx = 0. (B.24) 


For any function F{x), define 

oo 

[F] = J F{x)(/)'^{x)dx. (5.25) 



From (B.24), we have 

£n = . (i3.26) 

[Jn-l\ 

It is convenient to regard (B.21) as an electrostatic analog problem with — ^/^ 
as the electrostatic field, 

K = (l)'^ (5.27) 

as the dielectric constant, 

Dn = -\l^fn (5-28) 

the displacement field and 

(Jn={u- Sn)(f'fn-l (5.29) 
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the electrostatic charge density. The electrostatic field equation is 

D'n = On. (B.30) 

and (B.24) gives the condition for zero total charge, 



00 



y On{x)dx = 0. (B.31) 



Consequently, 

DniO) = 0. (B.32) 
At X = oo, because 0(oo) = 0, we also have 

Dnioo) = 0. (B.33) 

Furthernnore, Dnix) is positive and therefore 

ax) < 0. {BM) 

From (B.22), (B.28) and (B.30), we find 

X oo 

Dn{x) = j (7n{z)dz = - j an{z)dz (B.35) 

X 

and, on account of (B.23), 

oo oo 

Ux) = 1 - 2 / (t)-\y)dy j Gn{z)dz (5.36) 

X y 

which is equivalent to 

X y 
fn{x) = /n(0) - 2 / r\y)dy I an{z)dz. 



Because u{x) satisfies u{x) > and u'{x) < 0, the hierarchy theorem 
applies. Therefore 

Sn^l < Sn (5.37) 

^ <0 (5.38) 



fn-l 
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and 



f 

J n 
fn-1 



' < 0. (5.39) 



Define 

9nix) = fn{x) - fn-l{x)- (5.40) 

forn > 2 
and for n = 1, 

ei = ei- i. (5.42) 
For clarity, we will present our analysis in the fornn of several simple theorenns. 



Theorenn 1 



and 



(/^)'<0, (^)'<0, (5.43) 

g'n < (^-44) 



Proof From (B.40), 



f- = ^-l and ^ = ^-1. (5.45) 
Therefore, (B.38) and (B.39) lead to (B.43) and therefore 

9n y ^ Jn_ _ 9nfn-l ^ q 
'fn-l' fn-1 fn-1 

Because f^_i < and Qn, fn-i both positive, it follows then 

9n < < 0. (5.46) 

Jn-l 

Theorem 2 

en > 0. (5.47) 

Proof For n > 2, (B.37) gives e„ > 0. For n = 1, (B.26) and /o = 1 lead 
to 
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As we shall prove 



N = (^i-l)[l] = [n-i]>^e-t^. (5.49) 



By using (B.9), we obtain 
Introduce 

4 

As X varies fronn to oo, ^ follows a path V, starting fronn ^ = e~^^, increasing 
to ^ = 1 when x = 1, and then decreasing to ^ = 0. Thus, 

1 1 Q I 

[u-t-] = f (f)\u - ^)dx = [ , (^-52) 

Divide the positive x-axis into three sections: 

1. X from to 1; (B.53) 

4 9 

correspondingly ^ fronn e~39 to 1, and 5*0 from | to 0. 

2. X from 1 to VS; (B.54) 
^ from 1 to e~i^ and 5*0 from to |. 

3. a: from Vs to oo; (5.55) 
^ from e^s^ to and 5*0 from | to oo. 

Define 

Eq.(B.52) can be written as 

[u-\] = l{9) - ll{9) - lllig) (B.57) 
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where 

2 

3 



lig) = J Fix4So))e-'a'^dSo > 



Uig) = J Fix+iSo))e-^''°dSo > (B.58) 



and 

oo 

III(p) = / Fix+iSo))e-^'''dSo > 0, 



with X-{So) and x+{So) given below. 

In accordance with (B.6), when x varies from —2 to 2, 5*0 ranges from to 
|. in this range, for each 5*0 there are three real roots of x that satisfy 

x^-3x = 3So - 2. (B.59) 

Let 

cos^ = l-|5o. (B.60) 

With ^ = at 5*0 = and x = 1. For each 5*0 within and | designate these 
three roots as 

x+(So) = 2cosq;+(S'o), ^-(*S'o) = 2cosa_(S'o) 

and xo{Sq) = 2cosao{So), (B.61) 
with 

a+iSo) = -60° + ^ 

a_(5o)=60° + ^ (5.62) 

aoiSo) = 180° + ^. 

The x referred to in (B.53) and (B.54)-(B.55) are respectively the above 
x_{So) and x+{Sq) together with its analytical extension to x+ > 2. 



33 



When g = 0, by using 

Sq —— 1 X 

and (B.63) 
(3 + x){l-x^) = {l + x)(4-{l + xf^ , 
we have from (B.52) and (B.57) 



1 °9 I 



= 1(0) - 11(0) - III(O) = i (5.64) 

Take the expression for III(p) in (B.58). Since in its integrand e~^^'^° < e~^^, 
we have 

lll{g) < 6-5^11(0). (5.65) 

Next consider 

2 

I(^) - II(^) = / (f(x_(5o)) - F(x+(5o)))e-2^^°d5o (5.66) 



in which, since x+{So) > X-{So), 

F{x4So)) - F{x+{So)) > 0. (5.67) 
In addition, because in its integrand 

we find 

l{g) - ll{g) > e-l^ (l(0) - 11(0)) . (5.69) 
Combining with (B.64)-(B.65), we derive 

I(^) - II(^) - III(p) > e-t^(l(0) - 11(0) - III(O)) > h-l^ (5.70) 

which leads to (B.49) and ei > 0, and thereby completes the proof of 
Theorem 2. 
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Theorem 3 

' [1] 

and for n > 2 

_ [{U - en-l)9n-l\ ^ [{U - \)9n-l\ . „ 71 ^ 

[fn-l] [1] • ^ ^ 

Proof The first equation for n = 1 follows from (B.42). For n > 2, 



_ p c _ [(^ ^n-l)/n-l] / „ r,rf\ 

— c-n — c-n-l — FT ] • i^-'^j 

[jn-l\ 

In accordance with (B.40) and (B.26), 

fn-l = fn-2 + 9n-l (5.73) 

and 

[{u-Sn-l)fn-2]=0. (5.74) 

Thus, the equality in (B.71) follows. Since for all n, 

en>ei> i (5.75) 

we establish also the inequality in (B.71). 



Theorem 4 



[{u - l)gn] ^ 1 



<— (14p,(0) + 3GO, (5.76) 



[1] 8^ 
where 

G„ = maximum of — g'^{x). (5.77) 
Proof As in (B.52), we write 

O 1 oo 

= -7r9ni0)e-l' + ^ / <P\Agr, - Bg'Jdx (5.78) 
^g og ^ 



with 

{l + xy\l+x 3 + xJ {1+x) 
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and 

J^(^) = 777^ ^ 3. (i?.80) 

(1 + a;)"* 

Since the first term on the right hand side of (B.78) is negative, Theorenn 4 is 
proved. 

Theorenn 5 For < x < 1, 

-g'Sx) < ^ Pn-i(O). (5.81) 
Proof Integrating (B.21) fronn to x, and using (B.22), we obtain 

-U\^)f'Ax) = Jct^\z){u{z)-en)fn-l{z)dz 
^ 

and (5.82) 
1 ^ 

-^0^(^)/n-l(^) = / 0^(^)(^(^) - Sn-l)fn-2{z)dz. 



The difference between these two equations gives 

1 ^ 

-2^^(^)^n(^) = / 0^(^){(^(^) - Sn)gn-l{z) - enfn-2{z)^dz 





X 



< J (P\z)iuiz) - en)9n-liz)dz, (5.83) 


since -enfn-2{z) < 0. Fronn (B.44) and (B.49) it follows then 

9n-i{z)<gn-m, (5.84) 

Sn>Si>^ (5.85) 

and (B.83) becomes 

1^1 
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X ^ 

<9n-l{0)J(l)\z)(u{z)--yz 





X 







Because 
(B.86) leads to 



35^n-l(0) ^-2flgn(x) 



< ^ e 
2p 



By using (B.5), we derive 



-S'ni^) < ^i^a„-i(0) (B.87) 
which for < X < 1 leads to Theorenn 5. 

Lemma When n = 1 and x > 1, 

(-(l + a;)V;)'<0 (R88) 

provided 

g>2. {B.89) 

The prove of the lemma is given at the end of this Appendix. We will now 
proceed assuming its validity. 

Theorem 6 For x > 1 and g > 2, 

(-(1 + ^)V;)'<0 {B.90) 



and 



il + xYg'^) <0. (5.91) 
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Proof Since 



-(1 + ^ffL = 7^-4^---|(-(l + ^?n) (B.92) 

Jn-l Jn-2 Jl^ ^ 

and 



-(1 + ^fg'n = il^-§^---|(-(l + ^ffX (B-93) 

/n-1 /n-2 /l ^ ^ 

Using (B.39), (B.43) and the lennnna, we see that the derivatives of (B.92) and 
(B.93) are negative. 

Theorem 7 For g>2 and x>\, 

Gn = max ( - < - ^„_i(0). (B.94) 

Proof For x > 1, 

( - < (B.95) 
on account of (B.91); therefore 

Together with Theorem 5, (B.94) is proved. 
Theorem 8 For g > 2 and for all x > 0, 

9 /9\n 
9nix) < QniO) < - Qn-M <■■■<(-) {3.97) 



9 ^9 
Proof Since 

1 

9niO) = 9n{l) - J 9ni^)dx < + (B.98) 



and 

oo oo 

9nil) = - / 9ni^)dx = I [-(1 + ^)'9n]^^^^2 ' (^-99) 

On account of (B.91), for x > 1 

-(l + x)2^;(x)<-4^;(l); (5.100) 
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therefore (B.99) yields 

dx 

Combining this result with (B.98), we derive 



9n{l) < -4^;(1) / 7T--T2 = -2^Ul) < 2^,. (5.101) 



^„(0) <3Gn<^ ^n-l(O), (5.102) 



on account of (B.94). 
Theorem 9 For g > 2, 
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Proof Using (B.76), (B.94) and (B.97), we derive 



5 /9\n 

en < 7T7 ( - ) (^-103) 



^^^^^|^<^(14^n-i(0) + 3G„_i) 

1 / /9\«-i /9\"-i\ 15 /9\n-i 
<— 14- + 



8g\ \gy \g^ / 8g\g 
5 /9\n 
~ 24 V 

Substituting this result into (B.71), we complete the proof of Theorem 9. 

We will now turn to the proof of the lemma. For n = 1, (B.29) and (B.36) 
lead to 

oo 

f[{x) = 2(f)-\x) J {u{z) - ei)(j)'^{z)dz. (5.104) 

X 

Define 



F{x)^-]^{l + xff[{x) (5.105) 

^{x) ^ jP^, = 7r^e-2^^o(^) (5.106) 
^ ^ 4(l + a:)2 (l + a;)4 ^ ^ 



and 



rj{x) = / (u(z) - ei)cl)\z)dz. (5.107) 



X 
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We find 
with 



l' = 




r = 




p = 


{l + x 




4 


Q = 


l + x 



and q = ,^ + 2g{x'-l). 

Since > ;^ in accordance with (B.49), for a; > 1, we have 



u{x) = ^ < 81. 



1 

{x + iy 

Therefore 

r]{x) > 0; 

in addition, 

p{x) > and q{x) > 0. 
Regarding 77 = r]{^), we define 

Its derivative is 

dL (Pij pq' — qp' 

'd^^^de^ q' ■ 
Likewise, by differentiating (B.108), we derive 

At a; = 00, (B.106) and (B.107) yield 

^(00) = and 77(00) = 0; 

their ratio is 

^(00) q(oo) 2^^^- 
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From (B.115), we see that at a: = oo 

= m = 0. 

V / x—oo 



(5.120) 



As X decreases, since ^' and t]' both are negative, ^{x) and r]{x) increase and 
are both positive. Next, we shall prove that for a: > 1 and g > 2 



.e.. 



1 

Kx) = -^{pq' -qp') > 0; 



dL 4 , ^ 



L>0 



and therefore, because of (B.117), 



dF 



> 0. 



In turn, because of (B.llO), for a: > 1 

dF 

and that leads to, on account of (B.105) 

which is the lemma (B.88)-(B.89). 

To establish (B.121), we use the identity 



oo 



02 (x) 



(1 + x)" 

-1 oo 



— -^-\dx = ^ 



+ 



n+ 1 
22(1 + 



2"^ 

/? + 3 



e 3i 



+ 



2--ip 

11 + 2 



+ ••• + 



(l + x)"+2 2(l+x)"+i 
4 



)n-l 



(1 + x)^ 



(B.121) 

(B.122) 
(5.123) 

(5.124) 



(5.125) 



(5.126) 
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It follows then 

1 1\ , 3 



/ ^ - ^)dx = -—e 3^ 



(l + x)2 4/ 2g 
1 5 4 





00 



4^^ ' 'V(l + x)4 2(l + a;)3 
°°^2/ 1 1 9 \^ /3 103 







and 



where 



_J_7,2( 54 45 36 X 

+ 2V / ^"^^ HI + x)6 + 2(1 + x f ^ 22(1 + a;)4 + 2^(1 + x)^/ 

(B.128) 

[11 ^ _ 3 103 2403\ 
1 T ,o, ,/ 315 712 938 1020 



2Vi 'Hl + a;)8 2{l + xy 22(l + a;)6 23(l + a;)5 

+ u + r^]dx. (B.130) 

24(l + a:)4 25(l + a:)3y ^ ^ 

Next, using (B.48), we can write the above expression (B.129) as 

^' = i + i + J (^-131) 

where 

85 ^ , , 

7 = ^ + - (B.132) 

with 6 defined by (B.130). From (B.111)-(B.112), we find that X{x) defined 
by (B.121) is given by 

X(x) = gSi(x+lf-gx + —^-—^-3ei. (5.133) 

(1 + X)''' 
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Its derivatives are 

X'ix) = 2geiix + l)-g- — ^ (B.134) 

and 

X"ix) = 2g8i + ^^-^ > 0. (B.135) 

At X = 1, 

>:{l) = Ag8,-g-\. (5.136) 

Using (B.131), we find 

y(i) = A + ^. (5.137) 

Neglecting 0(e~3^), (5 is positive; therefore 7 > and A'(l) > 0. Because 
\"{x) > 0, \'{x) > for all x > 1, and the nnininnunn of \{x) is at x = 1 with 

In order that A(l) > 0, we require 

647 > g[- ^ . (5.139) 



4^-3 

Assuming g not too small so that we can neglect 5/g in (B.132); hence 

7 = f (5-140) 

and 

647 ^ — = 10.625. (5.141) 
8 

The righthand side of (B.139) is 

fl(^^) = f when 3 = 1.5 (B.142) 

consistent with the inequality (B.139). If in (B.132), we take 5/g into account, 
then a sufficient condition for A(l) > and therefore X{x) > for x > 1 is 

g>2. (5.143) 

This, together with (B.122)-(B.125) complete the proof of the lemma. 
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Table 1 



g = 3 








2 

9 max 


2 

9min 


7 

k 


a 






(- + 1)(- + - + 1) 


(! + l) 


0.5 


0.4 


.46 


.25 


10.69 


5.06 


1 


0.6 


.72 


.5 


11.56 


-7-10 

7.13 


1.5 


0.8 


.98 


.75 


11.86 


8.29 


2 


1.2 


1.23 


1.0 


9.33 


7.13 


2.5 


1.3 


1.49 


1.25 


10.79 


8.52 


3 


1.6 


1.74 


1.5 


10.06 


8.29 


3.5 


1.8 


1.97 


1.75 


10.31 


8.64 


4 


2.1 


2.24 


2.0 


9.82 


8.41 



Figure caption 

Fig. 1. Iterative radial wave functions ^n(^) and their corresponding en- 
ergies En for both the boundary conditions /n(0) = 1 (upper curves) and 
fn{oo) = 1 (lower curves). The parameters are g = 3, k = 2 and a = 1.2. 

Fig. 2. Final 7^(r) = 7^iv,/(0 and = E^v,/ that satisfy (2.11) for g = 3, 
/ = and = 3, 4, 5, 6. The overall normalization factor for these curves are 
determined by requiring f{oo) = 1.1, 1.1 , 1.0 and 0.7 for the corresponding 
iV = 3, 4, 5 and 6. 

Fig. 3. Final 'TZ{r) = ^n,i{'>^) and E = E^j that satisfy (2.11) for g = 3, 
N = 3 and / = 0, 1, 2, 3. The overall normalization factor for these curves 
are all determined by requiring f{oo) = 1. 

Table caption 

Table 1. List of parameter a used in Figures 1-3 for different k = I + ^{N — 
1), but with the same (/ = 3. The two last columns are = (^+l)(^+^+l) 
and g^^^ — ia~^ different pairs (A;, a) used in these figures. 
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